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Abstract 

The electric losses in a bulk or film superconductor ex- 
posed to a parallel radio-frequency magnetic field may have 
three origins: In homogeneous vortex-free superconductors 
losses proportional to the frequency squared originate from 
the oscillating normal-conducting component of the charge 
carriers which is always present at temperatures T > 0. 
With increasing field amplitude the induced supercurrents 
approach the depairing current at which superconductivity 
breaks down. And finally, if magnetic vortices can pene- 
trate the superconductor they typically cause large losses 
since they move driven by the AC supercurrent. 

INTRODUCTION 

The phenomenon of superconductivity was discovered in 
1911 by Heike Kamerlingh-Onnes in Leiden. After he had 
achieved to liquify helium at the temperature of T = A.2K 
he observed that the resistivity of Hg became unmeasur- 
ably small below some "critical temperature" T c = 4.15 
K. A sensitive method to measure the residual resistivity 
in this "superconducting state" is to observe the temporal 
decay of the persistent "supercurrents" in a ring, say of Pb 
(T c = 7.2 K), Sn (T c = 3.72 K), or Nb (T c = 7.2 K) 
by monitoring the magnetic field generated by the circulat- 
ing current. It turned out [ 1 1 that the supercurrent does not 
decay measurably, even after several years. Ideally loss- 
free superconducting wires may thus be used to build coils 
which keep their magnetic field for years, after the wind- 
ings have been loaded with current and then are cut short 
by a superconducting switch. 

Thus, DC currents in a superconductor can flow practi- 
cally loss-free if they are not too large. However, it turned 
out that alternating currents (AC) in superconductors are 
not completely loss-free, in particular at high frequencies 
(RF = radio frequencies, MW = microwave frequencies). 
There are essentially three effects which cause energy dis- 
sipation during current flow in superconductors: 

(a) Even in ideally homogeneous bulk superconductors 
an electric field E cx uj (with frequency u/2ir) is required 
to accelerate the "superconducting electrons", the Cooper 
pairs of the microscopic BCS theory J2|. This electric field 
also moves the "normalconducting" electrons that are al- 
ways present at finite temperatures T > 0. The dissipated 
power of this effect is cx E 2 cx ui 2 . 

(b) When the current density inside the superconductor 
reaches the depairing current density jdp, the supercon- 
ducting order parameter is suppressed to zero at the place 
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where j — jd p . This means superconductivity disappears 
and electric losses appear. This nucleation of the normal 
state typically occurs at the specimen surface or in the cen- 
ter of Abrikosov vortices. In particular, when an increasing 
magnetic field H a is applied along a superconducting half 
space x > one initially has j(x) = (H a /X) exp(— x/X) 
when j <C jdp, but when H a reaches the thermodynamic 
critical field H c = $o/(v / 87rA^/io) one has j ps jd p near 
the surface, thus one has jdp ~ H c /\. Here <t>o = h/2e — 
2.07 • 10~ 15 Tm 2 is the quantum of magnetic flux, A is 
the magnetic penetration depth, and £ is the superconduct- 
ing coherence length. Within Ginzburg-Landau (GL) the- 
ory, valid near T = T c , the lengths A(T) = <(T) cx 
(T c - T)- 1 / 2 diverge as T -> T c , but the GL parameter 
k = A/£ is independent of temperature T. 

(c) Large dissipation may be caused by vortices inside 
the superconductor. These move under the action of the in- 
duced AC current, which exerts a Lorentz force on the vor- 
tices and causes them to oscillate and dissipate energy. At 
high frequencies the amplitude of this oscillation is smaller 
then the range of possible pinning forces caused by ma- 
terial inhomogeneities, e.g., precipitates or defects in the 
crystal lattice. This vortex dissipation then cannot be sup- 
pressed by introducing pins. 

For a flat bulk type-II superconductor (defined by k > 
l/v2) in thermodynamic equilibrium it is favorable that 
part of the magnetic flux penetrates in form of Abrikosov 
vortices when the applied field H a equals or exceeds the 
lower critical field H cl sa ($ /47rA 2 /xo)(ln/c + 0.5) (for 
K > 1.5) and has not yet reached the upper critical field 
H C 2 pb $o/(27r£ 2 /io) where superconductivity vanishes. 
The penetration of vortices at an ideally flat surface may 
be delayed by a surface barrier leading to a higher pene- 
tration field Hp pa H c > H c \ (overheating). On the other 
hand, with superconductors of finite size, demagnetization 
effects may allow the vortices to penetrate already at much 
lower fields. In particular, for a large film of width w and 
thickness d<w, the penetration field is strongly reduced, 
H p /H c i ~ d/w . . . y/d/w <C 1 depending on the edge 
profile, see below. An infinitely large thin film will thus be 
penetrated by any perpendicular magnetic field component, 
even if very small. 

AC RESPONSE OF VORTEX-FREE 
SUPERCONDUCTORS 

The puzzling fact that superconductors may carry loss- 
free DC current but AC currents exhibit electric losses, was 
explained by the two-fluid model of Gorter and Casimir in 
1934 [1]. Later, the microscopic BCS Theory [2j essen- 



tially confirmed the two-fluid picture, giving for its phe- 
nomenological parameters a microscopic interpretation and 
explicit expressions. 

Two -Fluid Model 

The phenomenological two-fluid model assumes that the 
total electron density is composed of the density of super- 
conducting electrons n s and that of normal electrons n n , 
which have different relaxation times r s and r„. Histori- 
cally, Gorter and Casimir assumed n n cx t 4 (t — T/T c ) 
and n s cx 1 — i 4 . As in the Drude model JT|, the drift ve- 
locity v of each of these two fluids should obey a Newton 
law, 



m dv/dt — eE — mv/r 



(1) 



with to and e the mass and charge of the electron. The total 
current density J = J s + J n is the sum of the supercurrent 
J s = en s v s and the normal current J n = en„v„. (In all 
other Sections of this paper the current density is denoted 
by j.) If one assumes t s = oo one obtains from Eq. (1) the 
first London equation 



dJ s /dt = (n s e 2 /ro)E = E/(/z A 2 ) 



(2) 



with A = (m/n s e 2 /xo) 1 / 2 the London depth. In the Lon- 
don gauge where E = —dA/dt (induction law) Eq. (2) 
may be written in form of the second London equation 
Js = (MoA 2 )~ A. For the normal electrons one may as- 
sume r n <l/u with periodic electric field E cx exp(iu>t). 
This gives for the normal current 



(n„e 2 T„/m)E. 



(3) 



Defining the complex conductivity o(uS) = a\(uf) — £02(0;) 
by J = <j{uj)E cx exp(jwi), one obtains 



cti(cj) = {-Kn s e 2 /rauj) 5{oS) + n n e 2 T n /m, 
02(0;) = n s e 2 /mijj = (/loA 2 ^) -1 3> a±. 



(4) 



For a normal conductor this yields <j\ = a n , a 2 = 0, and 
the skin depth 5 s ia n = (2/ fioOnU)) 1 ' 2 . For superconduc- 
tors at u) = the <5-function in o\ reflects the ideal DC 
conductivity, while at finite frequencies the inductive part 
dominates, 02 ^> o\ . However, the small dissipative part 
ci = n n e 2 T n /m is important since it causes the AC losses. 
In the situation with an AC magnetic field parallel to the su- 
perconductor surface, the current is forced (a current bias 
as opposed to a voltage bias) and one has a dissipation per 
unit volume pj 2 = Rc{l/a}J 2 (cri/cr|)J 2 rj a x E 2 
since a\ <C 02. The dissipation is thus proportional to 
n n uj 2 . 

The sum a = o\ — 102 is analogous to a circuit of 
a resistive channel 1/R cx <j\ in parallel to an induc- 
tive channel of admittance 1/iuiL cx a%. Below a fre- 
quency ojq = R/L this circuit is mainly inductive and 
above mainly resistive. The ratio of the currents in the two 
channels is J s /J n = n s J (n n ujT). This defines a crossover 
frequency a; ps (n s /n„)(l/r„) w {n s /n n ) ■ 10 11 Hz (TJ. 



When the superconductor forms the inner wall of a mi- 
crowave cavity with incident parallel magnetic field of am- 
plitude H inc , this wave is almost ideally reflected by the 
wall since a surface screening current J s = 2H mc is in- 
duced. The small dissipated power per unit area is then 
P s = J 2 R S , where 



(TV/lcrl 2 « S^CTi/a 2 



(5) 



is the surface resistance (e.g. in units f2). Here 5 = 
[2/fio(\a\ + o'2)w] 1 / 2 is the general skin depth reducing 
to the superconducting penetration depth A or to the skin 
depth (5 s ki n in the super or normal conducting limits. For 
the superconducting wall one has R s w (7i/1qA 3 cl> 2 /2 and 
the absorbed versus incident power of this wall is [ 1 1 
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(6) 



The quality factor Q of the superconducting cavity is thus 
inversely proportional to R s cx Q^ 1 cx n n u> 2 . 

Microscopic Theory 

After the BCS theory [2 1 had given the microscopic ex- 
planation of superconductivity, the complex AC conductiv- 
ity was calculated within this weak-coupling theory [3] O . 
In the extreme local limit (A <C £0 = hvp/n/S. with vf 
= Fermi velocity and A = energy gap; this assumption ac- 
tually is not satisfied for type-II superconductors with GL 
parameter k > 0.7), in the impure limit (electron mean 
free path I = vpr <C £o)> an d f° r frequencies below the 
energy-gap frequency {hv < 2 A, v = u/2ir) the resulting 
AC conductivity may be expressed as two integrals over an 
energy variable, 



2£= / h, 2 (e,A,T,u)de. 

<Jn J 



(7) 



where a n = ne 2 T/m = ne 2 l/pp (pp — mvp = Fermi 
momentum, n = electron density) is the Drude conductiv- 
ity in the normal state and fx and f-2 are some functions. 
Evaluating these integrals for the case u <C T <C A (in 
units Ti — kg = 1) one obtains [5] 

o 1 _ 2A ex 7 AN h n 9T 
V n ~ ~ exp V t) n 4w ' 

^ = ^. (8) 

The dissipative part <j\ and inductive part a% may be written 
in the form of the two-fluid model: 



o\ 

<T 2 



n qp e 2 l/p F . 
n s e /mui , 



(9) 



where n qp is the quasiparticle density (replacing the nor- 
mal electron density n n of the two-fluid model) and n s the 
superconducting electron density, 
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(10) 



The quality factor Q of the resonator is now 

Q" 1 oc R s R3 ^1\ 3 (Tiuj 2 oc n qp u 2 . (11) 

Since the quasiparticle density n qp ~ exp(— A/T) 
strongly decreases at low temperatures T, Q should in- 
crease drastically. Note that with increasing purity (in- 
creasing I) u\ increases but the penetration depth A w 
Apurei/l + Co A decreases in Eq. (11). Thus, maximum 
Q is reached at some intermediate, not too high purity of 
the superconductor. 

High-Purity Niobium 

For the high-purity Nb used in the TESLA cavities, 
the frequency dependent surface resistance has been com- 
puted by Kurt Scharnberg within the Eliashberg model that 
extends the BCS model to strong coupling superconduc- 
tors 0. Strong (electron-phonon) coupling effects change 
the amplitude and the temperature dependence of the gap 
parameter, they lead to a renormalization (enhancement) 
of the quasiparticle mass, which in turn affects the Lon- 
don penetration depth, and they result in temperature and 
energy dependent quasiparticle lifetimes. The electron- 
phonon interaction enters in form of the Eliashberg func- 
tion a 2 F(uj) which was taken from tunneling experiments. 
A Coulomb pseudopotential /i* = 0.17 and a Coulomb 
cut-off u) c = 240 meV were used. At sufficiently low T 
and low oj of the incident radiation, inelastic scattering is 
negligible and only disorder induced elastic scattering is 
important, which is parameterized by the normal state scat- 
tering rate Tn = l/2r. This is fit to the surface resistance 
R s measured at v = 1.3 GHz, yielding Tjv ~ 1 meV and 
t ~ 3 • 10~ 13 sec. Nonlocal effects (wave vector q > 0) 
were disregarded, which is partly corrected for by using the 
lifetime r fitted at 1.3 GHz. 

With these assumptions the surface resistance R s m 
^cti/1qA 3 w 2 of high-purity Nb was computed at T = 2 K. 
Note that R s is related to the reflectivity r of the metal by 
R s = (Z /4)(l - r) where Z = (eoc) -1 = 377 Q is 
the impedance of the vacuum. Starting from R s ps 20 nfi 
at 1.3 GHz the resistance rises to a few /ifl at 600 GHz 
and then exhibits a large step at 750 GHz to a value of 15 
mfL Above this energy-gap frequency v = 2A/h one has 
nearly constant R s till at least 2000 GHz. 

VORTICES IN SUPERCONDUCTORS 

Ginzburg- Landau and London Theories 

Before the microscopic explanation of superconductiv- 
ity was given in 1957 by BCS [2 | there were very power- 
ful phenomenological theories that were able to describe 
the thermodynamic and electrodynamic behavior of super- 
conductors. In 1935 Fritz and Heinz London established 
the London theory, see Eq.(2) above, and in 1952 Vitalii 
Ginzburg and Lev Landau conceived the Ginzburg-Landau 
(GL) theory. The GL theory may be written as a variational 
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Figure 1: Magnetic field B(r) and order parameter | - 0( r )| 2 
of an isolated flux line calculated from Ginzburg-Landau 
theory for GL parameters k = 2, 5, and 20. Note that the 
field in the vortex center is B(0) s» 2B C \. 
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Figure 2: Two profiles of the magnetic field B(x, y) and 
order parameter \ ^f(x, y) | 2 taken along the x axis (a nearest 
neighbor direction) for vortex lattices with lattice spacings 
a = 4A (bold lines) and a = 2A (thin lines). The dashed 
line shows the magnetic field of an isolated flux line from 
Fig. 1. Calculated from GL theory for n = 5 [ 10 1. 

problem that minimizes the spatially averaged GL free en- 
ergy density, 

<- m + i|^| 4 + KzCV + A^I 2 + (AVx A) 2 ) 

= Minimum . (12) 

Here ip(r) is the complex GL-function, or order param- 
eter, and A(r) the vector potential of the magnetic in- 
duction B = V x A. The two lengths are the mag- 
netic penetration depth A (usually taken as unit length) 
and the GL coherence length £; both lengths diverge as 
the temperature T approaches the critical temperature T c , 




Figure 3: Current stream lines, coinciding with the con- 
tours of B(x, y) and \ip(x, y)\ 2 . Abrikosov solution for the 
ideal vortex lattice near the upper critical field B C 2- 




Figure 4: Vortex lattice made visible by decoration with 
iron micro-crystallites. Top: Nb disk, 1 mm thick, 4 mm 
diameter, T = 4 K, B a = 985 Gauss, vortex spacing 
a = 170 nm (U. Essmann and H. Trauble 1968). Bottom: 
YBa 2 Cu307_«5, T = 77 K, B a = 20 Gauss, a = 1200 nm 
(D. Bishop and P. Gammel 1987). 

A oc £ oc (T c - TY 1 ! 2 . Their ratio, the GL parameter 
k = A/£ within GL theory (valid near T c ) is independent 
of T. 

The GL theory can be derived from the microscopic BCS 
theory (L. P. Gor'kov 1959) in the limit T c - T < T c , 
yielding for the GL function ip(r) = A(r)/ABcs where 
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Figure 5: Ideal reversible magnetization curves of a long 
superconducting cylinder or slab in parallel field B a com- 
puted from GL theory iflOll . 
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Figure 6: Magnetic field lines of vortex lines in and near 
a superconductor of finite size, and the circulating super 
currents (schematic). 

A is the energy gap function. The London theory follows 
from GL theory in the cases when the magnitude of the 
order parameter is nearly constant, \tp\ 1. This condi- 
tion is fulfilled when £ is small as compared to the spec- 
imen extension and to A, requiring k ^ 1. An arrange- 
ment of straight or arbitrarily curved vortex lines positioned 
at v v (z) — [x v {z),y v {z),z] (v — 1, 2, 3 . . .) then has a 
magnetic field that obeys the London equation modified by 
adding 5 functions centered at the vortex cores, 

(1 - A 2 V 2 ) B(r) = $ [ dr » S ^ r - r ") • (13) 

Ideal Vortex Lattice 

In 1957 Alexei Abrikosov, a thesis student of Lev 
Landau in Moscow, obtained a periodic solution of the 
Ginzburg-Landau equations and recognized that this cor- 
responds to a lattice of vortices of supercurrent, circulat- 




Figure 7: Magnetic field lines of a single vortex in a super- 
conducting film of thickness d = 8A (or half space z < 0). 
Analytical solution of London theory ifTTIl . 
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Figure 8: Magnetic field lines of the periodic vortex lattice 
in films of thicknesses d = 4A, 2 A, A, and A/2. From GL 
theory for « = 1.4 and B/B c2 = 0.04 JT2). The dashed 
lines mark the film surfaces. x\ = vortex spacing. 



ing around each zero of the order parameter and carrying 
a quantum of magnetic flux $ ; these vortex lines (or flux 
lines, fluxons) are energetically favorable when the applied 
magnetic field is between a lower and a higher critical field, 
H c \ < H C 2 (see Introduction and below). This solution ex- 
ists in bulk superconductors with GL parameter n > 1/ v2, 
called type-II superconductors. For this theoretical discov- 
ery Abrikosov received the Nobel Prize in Physics in 2003. 

Figure 1 shows the magnetic field B(r) and the order pa- 
rameter IV'tr)! 2 of one isolated vortex line for three values 
of the GL parameter k = 2, 5, 20. One can see that B de- 
cays over the length A and the vortex core has a radius s» £. 



For such not too small values of k to a good approximation 
the vortex field is the London solution, with the central sin- 
gularity smoothened over the core radius r c « \/2 £ 00, 



B v (r) 
K (x) 



K ( 



2ttA 2 " uv A 
ln(1.123/af), 



x < 1 



Y / 7r/2x exp(— x), x 3> 1 



(14) 



Kq(x) is a modified Bessel function. The interaction en- 
ergy of two vortices at a distance x 3> £ is f/mt = 
$oB v (x)/hq. 

Figure 2 shows cross sections of B(x, y) and \if)(x, y)\ 2 
along the nearest neighbor direction y = of the ideal tri- 
angular vortex lattice for two values of the average induc- 
tion B = (B) (with vortex spacing a = 2A and a = 4A) 
and k = 5. The dashed line is -B(r) for the isolated vor- 
tex. Figure 3 shows the contour lines of B(x, y) near B c2 
for the triangular vortex lattice. These lines coincide with 
the contours of \ip(x, y)\ 2 and with the stream lines of the 
supercurrents. 

The vortex lattice was first observed in the electron mi- 
croscope by U. Essman and H. Trauble [9| at our Max 
Planck Institute in Stuttgart, by decoration of the surface of 
a Nb disk with "magnetic smoke" generated by evaporating 
an iron wire in a He atmosphere of 1 Torr, see Fig. 4. The 
magnetization —A/ = B a — Softhe superconductor calcu- 
lated numerically ifTUl as function of the applied magnetic 
field B a is depicted in Fig. 5 for ideal (pin-free) long super- 
conductor cylinders or slabs in parallel B a (i.e., in absence 
of demagetization effects) with various k = 0.707 ... 10. 
At k = l/y/2 one has B c \ = B c = B c2 and the curve 
M(B a ) is the same as for type-I superconductors (with 
k < 1/V2), namely, B = 0, — M = B a (no penetrated 
flux) for B a < B c2 and B = B a , M = (complete pene- 
tration of flux) for B a > B C 2. 

When the superconductor is not a long cylinder or slab in 
parallel field, demagnetization effects shear the magnetiza- 
tion curves of Fig. 5 and reduce the field of first vortex pen- 
etration, see below. The vortices end at the upper and lower 
surface of finite-size specimens and send their magnetic 
field lines into the surrounding vacuum, see Fig. 6. The re- 
sulting modulation of B(x, y) just outside the surface can 
be observed by decoration (Fig. 4) and by magneto-optics 
or Hall probes. The magnetic field lines of one vortex in 
a thick film in a perpendicular magnetic field are depicted 
in Fig. 7 as obtained from London theory in [ 1 1 1. Figure 8 
shows the field lines of the periodic vortex lattice in films 
of thicknesses d = 4A, 2 A, A, and A/2 as calculated in lfl2l . 

Losses by Moving Vortices 

When a supercurrent flows in a superconductor, either 
applied by contacts or caused by a gradient or curvature 
of the local magnetic field, this current density j exerts a 
Lorentz force f = j x z$o on a vortex. The Lorentz force 
density on a vortex lattice is F = j x B. Neglecting a small 
Hall effect, the vortices move along this force with velocity 




moving vortex core induces electric field 
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Figure 9: Visualization of the origin of energy dissipation 
when a vortex moves with velocity v. Top: the dipolar elec- 
tric field lines induced by this motion run also through the 
normal conducting core (Bardeen-Stephen model). Bot- 
tom: During motion of the vortex core the order parameter 
relaxes (Tinkham term). 

v = r]~ 1 F where 77 is a drag coefficient or viscosity. The 
vortex motion induces an average electric field 

E = B x v = 7r X B x (j x B) = p ff j , (15) 
p s « {B/B c2 )p n . (16) 

Here pg is the flux-flow velocity, which at large average 
inductions B is comparable to the normal resistivity of the 
superconductor at that temperature (measurable by apply- 
ing a large field B a > B C 2). However, when only a few 
vortices have penetrated (B <§C B c i) one has much smaller 
resistivity pg -C p n . But even then the vortex-caused dissi- 
pation at low T is typically much larger than the dissipation 
caused by the normal excitations. 

Where does this resistive dissipation come from? There 
are two effects of comparable size, see Fig. 9. First, as 
pointed out by Bardeen and Stephen [HE), tne m °ti° n of 
the magnetic field induces a dipolar electric field that drives 
current through the superconductor and through the vortex 
core. If the vortex core is modelled as a normal conducting 
tube of radius r c s» £, the normal currents inside the vor- 
tex core dissipate energy that leads to the pg of Eq. (16). 
Second, as stated by Tinkham ifTl [T4l . the moving vortex 
core means that at a given position the order parameter \ tp\ 2 
goes down and up again when the core passes. If one as- 
sumes a delay of the recovery of \ip\ 2 by a relaxation time 
r « U/A one obtains an additional dissipation of the or- 
der of Eq. (16). These two sources of losses are nice for 
physical understanding. In the exact calculation of the dis- 
sipation of a moving vortex lattice from time-dependent GL 
theory [15] these two sources cannot be separated but the 
approximate Eq. (16) is essentially confirmed [ 16 1, also by 
microscopic theory 1 17]. The numerical and also the mea- 
sured flux-flow resistivity in the middle between the exact 



values and p n is somewhat smaller than the Eq. (16), i.e., 
for constant current source the real dissipation is lower. 
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Figure 10: Irreversible magnetization curves M(H a ) for 
a disk with radius a and thickness 26, b/a — 0.125, 
for various pinning strengths measured by the parameter 
aj c /H c i. The large hysteresis loop belongs to strong pin- 
ning aj c / H c i = 4. The small central loop is for the pin- 
free disk, whos vortex distribution is shown in Fig. 11. The 
reversible magnetization curve of a pin-free ellipsoid with 
same initial slope as M(H a ) is shown as dashed line. 



Pinning of Vortices 

When the material is inhomogeneous on the microscopic 
length scale of the vortex core £, then the vortices are 
pinned and cannot move as long as the Lorentz force does 
not exceed the pinning force, or the current density j is 
smaller than the critical current density j c , see the reviews 
lfl8l [191 l20l . In this way the electric losses caused by 
flux flow can be avoided, and completely loss-free conduc- 
tors of DC current can be tailored by introducing appropri- 
ate pinning centers into the material, e.g., precipitates and 
crystal lattice defects. For AC currents small losses remain, 
however. One source of AC dissipation is due to the (albeit 
small) concentration of normal carriers or excitations and 
can be understood from the two-fluid model as discussed 
above. The other source is the oscillation of vortices in the 
pinning wells. At small displacements u from their equi- 
librium position one may assume linear elastic binding of 
the vortices to the pins, with a force density — ku. Adding 
to this the viscose drag force — r\u and the Lorentz force 
one obtains the force balance equation in an AC current 
j QC oc exp(wji), 

j ac x B = -ku - iur]u . (17) 

One can see that at frequencies above fe/77, of order 
uj/2ir > 10 7 Hz, the viscose force dominates [21]. Pin- 
ning thus cannot prevent vortex oscillations at high fre- 
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Figure 1 1 : Penetration of vortex lines into pin-free cylin- 
ders with radius a and height 2b. Top: b/a — 2. Bottom: 
b/a = 0.3. Forced by the applied field H a , the vortices 
enter from the corners, but only when the applied field has 
reached some threshold field do they jump to the middle 
leaving a vortex-free zone near the surface. With further 
increasing H a the vortices eventually fill the cylinder uni- 
formly from the middle. This delayed penetration without 
pinning is called geometrical barrier. Such a barrier is ab- 
sent only for ellipsoid-shaped specimens. 

quencies. This vortex-caused dissipation increases as ui 2 , 
like the quasiparticle dissipation, cf. Eq. (11). 

Interesting theoretical problems are the statistical sum- 
mation of random pinning forces to obtain the critical cur- 
rent density j c , and the problem of thermally activated de- 
pinning |[T8l [T9 20 1 . The latter leads to vortex motion even 
at small currents densities j < j c due to finite temperature. 
This flux creep may be described by a highly nonlinear re- 
sistivity. In particular, a logarithmic dispersive activation 
energy for depinning, U(j) = Uo ln(j c /j), leads to an of- 
ten observed power-law current-voltage curve, 



E(j) = E eM-U{j)/k B T] = EoU/jc 



(18) 



with a creep exponent n — Uo/ksT. For n — 1 one has 
Ohmic behavior [free flux flow, Eq. (15)], for n ^> 1 one 
has flux creep, and in the limit n — > oo this power law 
yields the Bean model, in which j is either or j c : When 
at some position one has j > j c , the vortices rearrange im- 
mediately such that j is reduced to j c again. This concept 




Figure 12: Bean model with constant critical current den- 
sity j c for a superconducting bar with rectangular cross sec- 
tion 2ax2b{b/a = 0.35) put into a perpendicular magnetic 
field H a that first increases from to 0.5 (left column) and 
then decreases again (right column). The parameter 0.01, 
0.1, 0.3, 0.5, 0.3, 0.1, 0, -0.1 is H a /{aj c ). Shown are the 
magnetic field lines (solid lines) and the penetrating fronts 
(dashed lines) where the current density j (flowing along 
the bar) jumps from ±j c to (in the field-free and current- 
free core) or from j c to — j c (after full penetration of flux). 



is useful for DC currents and at not too high frequencies 
where the pinning forces exceed the viscose drag force. 

Geometry Effects 

The electromagnetic properties of a superconductor (and 
of any conductor or isolator) depend not only on the mate- 
rial but also on the geometry of the problem, i.e., on the 
shape of the specimen and on the way a magnetic or elec- 
tric field is applied. For example, the reversible magne- 
tization curves of a pin-free superconductor in Fig. 5 ap- 
ply to the unrealistic case of very long slabs or cylinders 
in exactly parallel field, where demagnetization effects are 
absent. For the still unrealistic situation of a perfect el- 
lipsoidal shape one may calculate from these ideal curves 
the reversible magnetization curves of any ellipsoid by us- 
ing the concept of the demagnetization factor. But when 




Figure 13: Bean model for the penetration of a perpendic- 
ular magnetic field B a = /JoH a into a thin rectangular film 
with thickness d <C width. The parameter H a /J c = 0, 
0.5, 1.5 measures H a in units of the critical sheet current 
Jc = djc- Shown are the stream lines of the sheet current 
in the film, 3(x, y) = J y, z) dz (left), and the contour 
lines of the magnetic field B z (x, y) in the plane z = of 
the film (right). 

the specimen shape is not an ellipsoid, then even for a pin- 
free superconductor the magnetization curves have to be 
computed numerically, since now the induction (or vortex 
density) inside the specimen is no longer spatially constant. 

It turns out that even without pinning such magnetiza- 
tion curves in general show a hysteresis, i.e., they are irre- 
versible and depend on the magnetic history, see Fig. 10. 
This irreversibility is due to a geometric barrier [22, 23] for 
the penetration of vortices as illustrated in Fig. 1 1 for cylin- 
ders (or long bars) with rectangular cross section: When the 
applied uniform field H a is increased, vortex lines enter at 
the corners, pulled by the screening currents (Meissner cur- 
rents) that flow at the surface, and held back by their line 
tension (like a rubber band). With increasing H a the vor- 
tices penetrate deeper and become longer. When the vor- 
tices from two corners meet at the equator, they connect 
and form one long vortex line that contracts and immedi- 
ately jumps to the specimen center. During this rapid jump 
all their elastic energy is dissipated by the viscose drag 
force F = r/v, see text above Eq. (15). With further in- 
creasing H a more vortices jump to the center, crossing the 
flux-free zone near the surface, and eventually the entire 
specimen is filled with vortices coming from the growing 
central zone. 

Flux penetration thus occurs with a threshold, over a "ge- 
ometrical barrier". The sudden onset of flux penetration to 
the center leads to the sharp maximum in the small (inner, 



pin-free) hysteresis loop of M (H a ) in Fig. 10. When H a is 
decreased again, the vortices leave the specimen essentially 
without barrier, and at H a = all vortices have left, i.e., 
one has B = and also M = (since no screening cur- 
rents flow anymore). The perpendicular field at which the 
first vortices enter at the corners of a pin-free long strip and 
a circular disk, both with rectangular cross section 2a x 2b, 
was computed in l23l : 

H pen P ~ Ha tanh ^0.36 b/a , 

H*lf w H c i tanh ^0.67 b/a. (19) 

In the presence of pinning the hysteresis loops of 
M(H a ) in Fig. 10 become larger. The area of such loops 
is the energy dissipated during one cycle due to depinning 
of vortices. When H c i is negligibly small as compared to 
H a , the hysteresis curves and the vortex density and cur- 
rents in a superconductor with pinning may be computed 
by treating it as a nonlinear conductor, Eq. (18). Figure 12 
shows how the magnetic field lines (and vortices) penetrate 
and exit a thick disk with pinning when an axial H a is first 
increased beyond the field of full penetration, and then is 
decreased again Il24ll25l . The chosen large creep exponent 
n = 50 practically reproduces the Bean model. 

Figures 10, 11, and 12 were computed by time- 
integration of an equation for the (scalar) current density j 
inside the superconductor; this method implicitly accounts 
for the infinitely extended magnetic stray field outside the 
specimen, without need to compute it and to cut it off. 
From the resulting current density the magnetic field lines 
are then easily calculated by the Biot-Savart law. 

A completely different geometry is shown in Fig. 13, 
namely, the current stream lines and the contours of the 
magnetic field B z (x, y) in a thin film or platelet of rect- 
angular shape 1 26 1 with pinning and large creep exponent 
n = 50 corresponding to the Bean model like in Fig. 12. 
An increasing magnetic field H a is applied perpendicular 
to the film. Initially, when H a <C J c = dj c is small, 
no magnetic flux penetrates the film, i.e., the circulating 
screening currents generate a magnetic field that in the film 
area is constant (of size —H a ) and exactly compensates the 
applied field H a . With increasing H a , magnetic flux pene- 
trates mainly from the middle of the sides of the rectangle 
(not from the corners), leaving still a flux-free zone in the 
middle. At and beyond some field of full penetration the 
current stream lines are concentric rectangles of constant 
distance, since the magnitude of the sheet current has sat- 
urated to the constant value J c = dj c . The magnetic field 
has then penetrated to the center, and the contour lines of 
B z (x,y) do not change anymore with further increasing 
H a . 

Penetration of First Vortex 

An important question for RF superconductivity is under 
what circumstances and at which applied magnetic field H p 
the first vortex enters the superconductor, since the pres- 
ence of even a few vortices can cause large losses. First 
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Figure 14: Gibbs free energy G of one vortex penetrating 
into a superconducting half space to a depth x, Eq. (22). Pa- 
rameter is the applied field H a in units of H&. The Bean- 
Livingston Barrier exists for H a < H c . 
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Figure 15: Nucleation of vortices as an arc of a circle at a 
planar surface (top), at a rectangular corner (middle), and 
at a rough surface (bottom, schematic). 



I summarize the expressions for the three critical fields 
which for type-II superconductors (with k > 1/ V%) obey 

B c i < B c < B C 2- 



B K = 
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4ttA 2 
$ 



(InK - 
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Figure 16: Supercurrents j z (x, y) in a bar with rectangular 
cross section 2a x 26 (b/a = 0.4) in the Meissner state 
with London penetration depth A = 0.025a. The currents 
(along the bar) are generated by a perpendicular applied 
uniform magnetic field H a \\z that penetrates to a depth A. 
Shown is one quarter of the cross section. Note the high 
(but finite) peak of j z (x, y) at the corners. The inset shows 
the magnetic field lines. 
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Figure 17: Field enhancement near the sharp edge of an 
ideal diamagnet. 
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While the thermodynamic (B c ) and upper (B C 2) critical 
fields are exact, the lower critical field B c i = ^qH c i has to 
be calculated numerically from the self energy of a vortex 
of length L, U sc \i — &qH c iL. The function a(re) is a good 
analytical fit to the numerical result of [ 10 1. At H a = H c i 
the nucleation of a vortex and motion to a depth x ^ A 
does not cost energy, see Fig. 14. However, the penetrat- 
ing vortex has to surmount a barrier such that the field of 
first penetration H p is larger than H c i- This barrier was 
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Figure 18: Enhancement of the magnetic field H at the 
equator of ideal diamagnetic cylinders, spheres, and strips 
of elliptic or rectangular cross section put in a uniform axial 
applied magnetic field H a . 

first predicted by Bean and Livingston (BL) [27| for a su- 
perconductor with planar surface in a parallel applied field 
H a . The Gibbs free energy G(x) for this case reads 



G(x) = $ (H a e 



±H v (2x) + (H cl -H a )). (21) 



-x/X 



In it the first term is the interaction of the vortex 
with the applied field H a or with its screening currents 
(H a / \)e~ x / x , the second term is the interaction with the 
image vortex (at position —x, of opposite orientation), and 
the third term is an integration constant. Using the fact that 
for not too small k one has B V (Q) w 2B C \ (see Fig. 1) one 
has with Eq. (14), B cX w (<S> /4ir\ 2 )K (r c /\) yielding 
with B c i (20) a core radius r c w £ exp[— «(«;)]. With this 
we may write G(x) (21) in the dimensionless form 



G(x) 



Ha_( p -x/\_ 



1 



-1- 



Ky^Ax^Wj/X) 
K (r c /X) 



(22) 



that is plotted in Fig. 14 for k ~ 1.3. Of course, this G(x) is 
a only approximate, in particular at small k, for which vor- 
tex penetration has to be computed numerically. Anyway, 
Fig. 14 shows that vortex penetration becomes favorable at 
H a = H c and that the Bean-Livingston barrier vanishes at 
H a K, H c > H c i. 

The assumption of BL that the entering vortex is long, 
straight, and exactly parallel to a planar surface is not very 
realistic. Alternatively, one may assume that the first vor- 
tex nucleates and penetrates in form of a small loop, say a 
half circle of radius R, see Fig. 15 top. The self-energy of 
this half circle is approximately ?7 se if = ttR^q/AttX 2 ^), 
putting the outer cut-off radius w R instead of A 3> R 
in the logarithm ln(A/£) — > ln(i?/£) « 1 when R is of 
order of £. The interaction of this vortex loop with the sur- 
face screening current of density j s is Uj B s» (7ri? 2 /2)$o j s 



(flux quantum times loop area times j s ). For a planar sur- 
face one has j s — H a /X directly at the surface. The crite- 
rion that Uj s > U se \f at H a > H p yields then 



(23) 



2-kXR R 

which is just the BL result. Thus, the assumption of a pen- 
etrating vortex loop does not change much the penetration 
field of a planar surface. 

However, when the surface has roughness with charac- 
teristic length > £, then vortices will penetrate at sharp 
points or cusps, see Fig. 15. At a corner with angle 
a = 90°, the screening current directly at the surface is 
strongly enhanced at this corner; Fig. 16 shows this for a 
superconducting bar with square cross section 2a x 2a and 
penetration depth A = 0.025a, to which a uniform trans- 
verse H a is applied. A rough estimate gives for this geom- 
etry an enhancement of the screening current at this corner, 
j s = CH a / A, by a factor C«4 The field of first vortex 
penetration H p is then reduced from Eq. (23) by just this 
factor, H p « H c /C w H c /4. 

For sharper corners the enhancement of j s and reduction 
of Hp are even larger. As shown in the textbook of Landau- 
Lifshitz (Electrodynamics of Continua) for an ideal dia- 
magnetic material at a corner with angle a (Fig. 17) the 
magnetic field diverges as H oc 1/r" with exponent f3 — 
(tt — a)/(27r — a), where r is the distance to the point 
of the corner. This gives H oc 1/r 1 / 3 for a — n/2 and 
H oc 1/r 1 / 2 for a -> 0. 

Similarly, an axially applied magnetic field flowing 
around an ideal diamagnetic cylinder, sphere, or disks 
with elliptical or rectangular cross section of aspect ratio 
b/a -C 1, is enhanced at its equator by factors 2, 3/2, a/b, 
or (a/6) 1 / 2 , respectively, due to the strong curvature of 
the field lines at this line, see Fig. 18. 

Vortices in thin films 

One has to distinguish two quite different types of vor- 
tices in thin film superconductors: vortices perpendicular 
or parallel to the film plane. In wide thin films with width 
w = 2a ^> thickness d = 26, the vortices will nearly al- 
ways run perpendicular across the film thickness, even in 
tilted applied field H a , because of the large demagnetiza- 
tion factor of this film. This means the circulating currents 
prefer to flow in the film plane. Only when H a is exactly 
parallel to the film surface, or when the film is coating a 
bulk superconductor that screens any perpendicular field 
component, then vortices parallel to the film plane may oc- 
cur. 

When the film is of finite size, one may use Eqs. (19) 
to estimate at which applied perpendicular field compo- 
nent H az the first vortices penetrate, namely already at a 
very small field, smaller than H c \ -J d/w. When the film 
edges are wedge-shaped or sharp, the penetration field is 
even smaller, cf. Fig. 17 and Fig. 18 (elliptical edge). Into 
infinitely extended or closed films (e.g., a Nb layer cover- 
ing the inner surface of a Cu cavity) any perpendicular field 



will penetrate since the field lines cannot flow around the 
film. Only when this film has holes or slits can some mag- 
netic flux cross the film via these holes, but the magnetic 
field in the holes will be larger than H az by at least the ra- 
tio of film area over the total area of all holes. However, 
the field in the holes will penetrate into the film when it is 
of the order of H c i times the square root of film thickness 
over hole distance. Thus, even such a perforated pin-free 
film will be penetrated by a perpendicular field that is very 
much smaller than H c \. The peaked magnetic field, cir- 
culating current, and pair interaction of perpendicular vor- 
tices in thin films were calculated for infinitely extended 
11281 and finite-size (e.g. rectangular) films [26 29, 30 1 . 

Pinning of vortices will not appreciably enhance all 
these penetration fields at high radio frequencies, where the 
(elastic) pinning forces are smaller than the viscose drag 
force. If the small applied perpendicular magnetic field is a 
DC field (e.g., some stray field or the earth magnetic field) 
then the additional RF field will even favor the penetration 
of the DC field in form of vortices, since it "shakes" the 
vortices. As shown in iFJTl l32l . shaking of vortices by an 
AC field oriented perpendicular to the vortices leads to the 
relaxation of irreversible currents if the AC amplitude ex- 
ceeds some threshold value. This vortex creep means that 
even in very small H az , perpendicular vortices will pen- 
etrate under the action of a large-amplitude RF field, and 
then these vortices oscillate and dissipate energy. 

The problem of parallel vortex lines in a thin film with 
d <C A was solved by Alexei Abrikosov (1964), Vadim 
Shmidt (1969), and in an elegant way by Alex Gurevich 
11331 . The lower critical field is enhanced in thin films as 
compared to bulk superconductors, 
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(24) 



and the field at which the surface barrier for vortex pene- 
tration disappears is also enhanced, 



R - $0 
p ~ 2^ 



(25) 



For example, a NbN film with £ = 5 nm, d = 20 nm 
has B c i — 4.2 T and B p — 6.37 T, much better than the 
penetration field B p rj B c = 0.18 T for Nb at low T. 

To enhance the operating RF amplitude in microwave 
cavities for accelerators and reduce the losses, Gurevich 
[34 1 suggests to use solid Nb or Pb with multilayer coating 
on its inner surface by alternating superconducting and in- 
sulating layers with d < A. This will prevent penetration of 
vortices into the bulk superconductor when the vortex pen- 
etration field Bp is large; e.g., for NbN films with d — 20 
nm the RF field can be as high as 4.2 T. From the elastic and 
viscose forces on a parallel vortex in a thin film, Gurevich 
estimates its characteristic relaxation time as 



2d^ A 2 /(£/dn 



(26) 



For a 30 nm Nb 3 Sn film this r ~ 10 12 s is much shorter 
than the RF period of 10~ 9 s. The maximum amplitude 



of the RF field at which the surface barrier of a single thin 
film coating disappears is of the order of the bulk H c of 
the film material, e.g., 0.54 T for Nb 3 Sn. Thus, Nb 3 Sn 
coating more than doubles the vortex penetration field for 
Nb, B p « B c = 0.18 T at low T. It appears that Nb 
cavities coated with a Nb 3 Sn layer or with NbN multilayers 
allow for much higher RF amplitudes than uncoated Nb, or 
Cu coated by a Nb film, if this can be achieved technically. 
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